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ABSTRACT 

Let  AT  be  a  3-D  convex  polyhedron  having  n  vertices.  A 
sequence  ^  of  edges  of  K  is  called  a  shortest-path  sequence  if 
there  exist  two  points  X,  Y  on  the  surface  S  of  K  such  that  ^  is 
the  sequence  of  edges  crossed  by  the  shortest  path  from  X  to  Y 
along  S.  We  show  that  the  number  of  shortest-path  sequences 
for  K  is  polynomial  in  n,  and  as  a  consequence  prove  that  the 
shortest  path  between  two  points  in  3-space  which  must  avoid 
the  interiors  of  a  fixed  number  of  disjoint  convex  polyhedral 
obstacles,  can  be  calculated  in  time  polynomial  in  the  total 
number  of  vertices  of  these  obstacles  (but  exponential  in  the 
number  of  obstacles). 


1.  Introduction 

In  this  paper  we  study  several  problems  related  to  the  problem  of 
calculating  the  Euclidean  shortest  path  between  two  points  in  3-dimensional 
space,  which  must  avoid  the  interiors  of  a  collection  of  polyhedral  obstacles 
having  altogether  n  vertices.  This  general  problem  seems  to  be  intractable, 
and  the  only  known  algorithms  for  it  require  exponential  time  ([SS],  [RS]), 
although  no  lower  bounds  are  known  as  yet  for  this  problem.  On  the  other 
extreme  hand  we  have  the  problem  of  finding  the  shortest  path  between  two 
points  in  3-space  which  must  avoid  the  interior  of  a  single  convex  polyhedral 
obstacle.   In  this  case  the  problem  is  solvable  in  time  0(n^\og  n)  ([SS],  [Mo]). 

Interpolating  between  these  two  extreme  cases,  one  might  consider  the 
problem  in  which  the  polyhedral  obstacles  consist  of  a  fixed  number  k  of 
disjoint  convex  polyhedra  (having  altogether  n  vertices),  and  attempt  to 
calibrate  the  complexity  of  this  problem  as  a  function  of  k  and  n. 

Baltsan  and  Sharir  [BS]  have  recently  studied  the  case  k  =  2,  and  have 
obtained  an  algorithm  for  this  case  running  in  polynomial  time  (almost 
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proportional  to  C>(n^log  n)).   However,  even  for  *  =  3  the  problem  becomes 

quite  difficult.   It  is  easy  to  see  that,  given  three  disjoint  convex  polyhedra  A, 

B,  C,  the  shortest  path  between  two  points  X,  Y  lying  respectively  on  the 

surfaces  of  the  polyhedra  A,  B,  which  also  passes  through  the  surface  of  C, 

must  consist  of  the  following  portions  (see  Fig.  1): 

(i)  a  shortest  path  along  the  surface  of  A  from  X  to  a  first  "take-off"  point  t^; 

(ii)  a  straight  segment  from  /i  to  a  first  "landing"  point  /i  on  C; 

(iii)  a  shortest  path  along  the  surface  of  C  from  Zi  to  a  second  take-off  point 

h; 

(iv)  a  straight  segment  from  t2  to  a  second  landing  point  /2  on  5; 

(v)  a  shortest  path  along  the  surface  of  B  from  I2  to  Y. 

It  is  also  easily  seen  that  the  take-off  and  landing  points  must  lie  on 
edges  of  the  corresponding  polyhedra  (unless  they  coincide  with  X  or  y 
themselves). 

If  one  knew  the  sequences  of  edges  on  A,  C,  and  B  crossed  by  the  three 
subpaths  (i),  (iii),  (v)  respectively,  then,  using  the  planar  unfolding 
procedure  described  in  [SS],  one  could  unfold  the  faces  of  A  (resp.  C,  B) 
crossed  by  the  subpath  (i)  (resp.  (iii),  (v))  so  that  they  all  lie  in  a  common 
plane,  and  so  that  the  corresponding  subpath  becomes  a  straight  segment. 
This  reduces  the  problem  to  that  of  finding  a  shortest  path  between  two 
points,  which  is  constrained  to  pass  through  four  specific  edges  (namely  the 
edges  containing  the  take-off  and  landing  points).  This  latter  problem  can 
then  be  solved  in  constant  time,  using  the  technique  described  in  [SS]. 

The  subpaths  (i)  and  (v)  start  or  terminate  at  fixed  known  points,  and 
the  results  of  [SS],  [Mo]  show  that,  for  a  polyhedron  with  n  vertices,  there 
are  at  most  0{n^)  sequences  of  edges  which  can  be  crossed  by  shortest  paths 
along  that  polyhedron,  emerging  from  a  fixed  source  point.  However,  the 
subpath  (iii)  starts  and  ends  at  unknown  points,  so  that,  a  priori,  there  could 
be  an  exponential  number  of  possible  sequences  of  edges  of  C  that  could  be 
crossed  by  the  subpath  (iii). 

In  this  paper  we  show  that  this  is  not  the  case,  and  that  there  are  only 
0{n^)  possible  sequences  of  edges  of  a  convex  polyhedron  K  with  n  vertices, 
which  could  be  crossed  by  shortest  paths  along  the  surface  of  K  between  any 
two  source  and  target  points.  This  implies  that  the  number  of  edge  sequences 
that  have  to  be  explored  in  searching  for  our  desired  shortest  path  is  only 
polynomial,  so  that  shortest  paths  amidst  three  convex  polyhedra  can  be 
calculated  in  time  polynomial  in  the  total  nimiber  of  their  vertices.  In  fact, 
this  observation  can  be  easily  generalized  to  yield  polynomial-time  algorithms 
for  finding  shortest  paths  amidst  any  fixed  number  of  convex  polyhedral 
obstacles. 

The  paper  is  organized  as  follows.  In  Section  2  we  prove  our  main  result 
about  a  polynomial  upper  bound  on  the  number  of  sequences  of  edges 
crossed  by  shortest  paths  along  a  convex  polyhedron.  In  Section  3  we  extend 
the  observations  made  above  to  obtain  polynomial  time  shortest-path 
algorithms  for  any  fixed  number  of  convex  polyhedral  obstacles. 


2.   On  the  number  of  shortest-path  edge  sequences  on  a  convex  polyhedron 

Let  /sT  be  a  3-D  convex  polyhedron  having  n  vertices.  Let  S  be  the 
surface  of  K.  For  each  pair  of  points  A,  5  in  5  we  denote  by  tt(A^)  the 
shortest  path  along  S  from  A  to  5,  if  this  path  is  unique,  and  by  ^{AyB)  the 
sequence  of  edges  of  K  crossed  by  tt(AJ5)  (if  A  or  5  lies  on  an  edge,  we 
indude  this  edge  in  ^(A,B)).  If  there  exists  more  than  one  shortest  path  from 
A  to  5,  7t(A,5)  will  denote  any  one  of  these  paths.  We  call  ^(A,B)  a 
shortest-path  edge  sequence  (induced  by  A  and  B).  Fig.  2  illustrates  these 
definitions. 

We  first  recall  a  few  basic  properties  of  shortest  paths  along  S  (for  more 
detail  of  which  see  [SS]).  Let  A,  B  be  two  points  on  S,  then  7t(A,5)  does  not 
pass  through  a  vertex  of  K  and  is  geodesic,  in  the  sense  that  it  is  a  polygonal 
path  which  bends  only  at  edges  of  K,  such  that,  with  each  edge  that  it  crosses 
it  subtends  equal  incoming  and  outgoing  angles.  B  is  called  a  ridge  point 
relative  to  A  if  there  exists  more  than  one  shortest  path  from  A  to  5.  A 
shortest  path  from  A  cannot  path  through  a  ridge  point  relative  to  A.  The 
collection  of  ridge  points  relative  to  A  consists  of  O(n^)  straight  segments 
linked  together  so  as  to  form  a  tree  whose  leaves  are  vertices  of  K,  and 
whose  nodes  of  degree  2  are  points  lying  on  edges  of  K.  The  ridge  points 
partition  each  edge  of  K  into  0{n)  intervals  such  that  for  each  such  interval  /, 
the  sequences  ^(A,C),  for  C  ^  /,  are  all  the  same. 

Lemma  2.1:  Let  A,  5,  C  be  three  distinct  points  on  S.  Then  either  one  of  the 
paths  Ti(AJB),  '7r(A,C)  is  a  subpath  of  the  other,  or  these  two  paths  meet  only 
at  A. 

Proof:  This  is  well  known  (cf.  [SS]),  and  is  proved  as  follows.  Suppose 
neither  of  the  paths  •7r(A,5)  and  '7r(A,C)  is  a  subpath  of  the  other,  and  that 
they  meet  at  some  point  D .  Then  the  lengths  of  their  portions  between  A  and 
D  must  be  equal,  so  that  D  must  be  a  ridge  point  relative  to  A,  contradicting 
the  property  that  shortest  paths  cannot  pass  through  a  ridge  point.  □ 

Lemma  2.2:  Let  A,  5,  C,  D  be  four  distinct  points  on  S.  Then  the  paths 
TT(ArB)  and  TriCJD)  intersect  in  at  most  one  point,  unless  either  one  of  these 
paths  is  a  subpath  of  the  other,  or  their  intersection  is  a  shortest  path 
between  one  of  the  points  A,  B  and  one  of  the  points  C,  D. 

Proof:  Suppose  that  the  special  cases  just  mentioned  do  not  occur,  and  that 
the  paths  'tt{AJB)  and  i7(C,£>)  intersect  in  two  points  X,  Y.  \\^thout  loss  of 
generality  suppose  that  Y  lies  on  it  (A  ,5)  between  X  and  B,  and  that  it  lies  on 
'ir(C,D)  between  X  and  D.  Then  the  shortest  paths  iriXyB)  and  'Tr(X,D)  are 
subpaths  of  the  paths  it  (A  ,5),  'it(C,Z))  respectively,  and  by  Lemma  2.1  they 
cannot  intersect  at  Y,  a  contradiction  which  proves  the  claim.  □ 

To  estimate  the  number  of  shortest-path  edge  sequences  along  S,  we  first 
consider  the  restricted  case  in  which  the  source  point  A  of  the  shortest  path 
lies  on  an  edge  ei  of  K,  and  the  target  point  B  lies  on  another  edge  e2-  Our 
approach  is  to  consider  the  2-D  space  Z  =  ei  x  62  of  all  possible  such  pairs 
of  source  and  target  points,  and  partition  it  into  regions,  such  that  for  each 
such  region  R  there  exists  a  small  set  C(R)  of  edge  sequences  such  that 
C(A,5)  ^  C{R)  for  all  (A ,5)  €  R. 


To  obtain  these  regions,  consider  first  a  fixed  pair  (A ,5)  €  Z  such  that 
neither  A  nor  B  is  a  ridge  point  relative  to  any  vertex  of  K.  For  each  vertex  X 
of  K  let  7(X)  denote  the  concatenation  of  (the  unique  paths)  irCA^)  and 
it(X,B).  Then  we  have 

Lemma  2.3:  (a)  Let  X  and  F  be  distinct  vertices  of  K.  Then  -/(X)  and  7(7) 
intersect  in  at  most  one  point  (in  addition  to  their  common  endpoints  A  and 
B),  unless  both  A  and  B  lie  on  Tr(X,y). 

(b)  For  each  vertex  X  of  /sT,  the  path  7r(A,5)  does  not  intersect  7(X)  (except 
at  the  points  A  and  B),  unless  either  A  lies  on  Tr(B^)  or  B  lies  on  7t(A,X). 
Proof:  (a)  Note  that,  by  Lemma  2.1,  the  paths  'it(A,X)  and  ir(A,y)  can  meet 
one  another  only  at  A,  because  'n{AX)  cannot  pass  through  the  vertex  Y,  and 
vice  versa.  Similarly,  the  paths  'rr(X,B)  and  7r(7,B)  can  meet  only  at  B.  By 
Lemma  2.2  the  paths  'it(A,X)  and  'n{Y,B)  intersect  in  at  most  one  point, 
unless  they  overlap  in  one  of  the  manners  stated  there.  Since  shortest  paths 
cannot  pass  through  vertices,  it  is  easy  to  check  that  the  only  way  in  which 
these  paths  can  overlap  is  when  both  A  and  B  lie  on  Tr(X,y).  Similarly,  unless 
this  special  case  occurs,  it(Z,B)  and  'n{A,Y)  also  intersect  in  at  most  one 
point,  so  that  altogether  ^(X)  and  7(7)  intersect  in  at  most  two  points. 

Even  though  this  result  suffices  for  subsequent  analysis,  we  can  in  fact 
show  that  7(X)  and  7(7)  cannot  intersect  in  more  than  a  single  point.  Indeed, 
suppose  that  'niAX)  and  'ir(B,y)  intersect  at  a  point  D,  and  consider  the 
closed  curve  8  =  'it(X,B)  II  'ir(B,Z>)  II  7r(D,X).  It  is  easily  checked  that  this 
curve  is  simple,  so  that  by  the  Jordan  curve  theorem  it  partitions  S  into  two 
disjoint  regions  R^,  /?2-  Furthermore,  the  two  paths  'Tr(A,£>)  and  Tr(y,D) 
cannot  intersect  8  (except  at  D),  and  since  they  are  portions  of  curves  which 
cross  one  another  at  D  (with  their  extensions  past  D  forming  part  of  8,  it 
follows  that  they  both  must  lie  on  the  same  side  of  8,  say  in  i?i  (see  Fig.  3). 
But  then  it  (A,  7)  cannot  intersect  'niXJB),  for  otherwise,  since  this 
intersection  must  be  transversal,  Tr(A,y)  would  have  to  intersect  8  at  another 
point,  which,  as  can  be  easily  checked,  is  impossible. 

(b)  Suppose  that  'ir(A^)  intersects  7(X).  Then  without  loss  of  generality  we 
can  assume  that  'Tr(A,B)  intersects  ir(A,X),  which,  by  Lemma  2.1,  is  possible 
only  if  B  lies  on  'ir(A,X)  {X  cannot  lie  on  'ir(A,B)  because  shortest  paths 
cannot  pass  through  vertices  oi  K).  o 

Let  us  now  draw  all  the  paths  7 (AT),  for  X  a  vertex  of  K.  They 
collectively  partition  S  into  0{rr)  (simply)  connected  regions,  and  Lemma 
2.3(b)  implies  that  'Tr(A,B)  is  wholly  contained  in  just  one  of  these  regions. 
Let  R  be  that  region;  the  boundary  of  R  then  contains  both  A  and  B,  and  we 
denote  by  81,  82  the  two  portions  of  this  boundary  connecting  A  to  B.  Note 
that  R  contains  no  vertex  of  K  (except  possibly  on  8^  and  82).  Hence  if  e  is 
an  edge  of  K  which  intersects  ir(A,B)  (that  is,  e  ^  ^{A,B))  then  e  must  also 
intersect  both  81  and  82,  where  each  of  these  intersections  can  also  be  at  an 
endpoint  of  e.  Conversely,  if  e  is  an  edge  of  K  which  intersects  both  81  and 
82,  then  it  is  clear  that  e  must  also  intersect  it(A,B),  i.e.  e  6  ^{A,B). 

Furthermore,  let  e  be  an  edge  which  intersects  both  81  and  82.  Let  e'  be 
a  subsegment  of  e  whose  endpoints  pi,  p2  ^^  ^^  Si»  ^2  respectively,  and 
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whose  interior  is  wholly  contained  in  R.  We  claim  that  if  7r(A,B)  is  indeed 
contained  in  R  then  it  must  intersect  e' .  This  is  because  e'  splits  R  into  two 
disjoint  subregions,  one  of  which  contains  A  on  its  boundary  and  the  other 
contains  B  on  its  boundary.  Thus  any  continuous  arc  connecting  A  to  B  within 
R  must  cross  e' . 

This  implies  that  if  'tt(A^)  is  contained  in  R,  then  the  subsegment  e' 
defined  above  must  be  unique,  for  if  e  contained  two  such  subsegments, 
'n{A,B)  would  have  to  intersect  both  of  them,  contradicting  the  fact  that  a 
shortest  path  cannot  cross  an  edge  of  K  more  than  once. 

Finally,  suppose  that  tt{AJB)  is  contained  in  R,  and  let  J,  e  be  two  edges 
of  K  intersecting  both  81,  82.  Let  d'  =  P1P2.  ^'  =  ^1^2.  be  their  unique 
subsegments  connecting  8^  to  82  as  defined  above,  withpi,  q-^  6  8^,  andp»2. 
^2  ^  82-  Then  we  claim  that  pi  is  nearer  to  A  along  81  than  q-^  if  and  only  if 
P2  is  nearer  to  A  along  82  than  ^2-  This  property  follows  by  general 
topological  considerations  from  the  fact  that  d'  and  e'  do  not  intersect  one 
another,  and  that  each  of  them  partitions  R  into  two  disjoint  portions,  so  that 
the  other  subsegment  is  contained  in  just  one  of  them.  Moreover,  similar 
considerations  imply  that  pi  is  nearer  to  A  along  8^  than  q^  if  and  only  if  d 
precedes  e  in  the  sequence  ^(A,5).  See  Fig.  4  for  an  illustration  of  these 
observations. 

All  these  observations  lead  to  the  following  procedure  for  calculating  a 
small  collection  C{A,B)  of  potential  edge  sequences,  which  ^{A,B)  can 
assume. 

PROCEDURE  EDGEJSEQUENCES  (A,B): 

(I)  Calculate  all  shortest  paths  -niAyX),  'n{XJB),  for  X  a.  vertex  of  K,  to 
obtain  the  curves  y{X).  Find  all  the  intersection  points  between  these 
curves,  and  obtain  the  decomposition  of  5  into  open  connected  regions 
by  these  curves. 

(n)  For  each  such  region  R  whose  boundary  contains  both  A  and  B ,  find  all 
intersections  of  edges  of  K  with  each  of  the  half -boundaries  81,  82  of  R 
connecting  A  to  fi . 

(HQFind  the  collection  S  of  all  edges  of  K  intersecting  both  81  and  82.  For 
each  edge  e  in  S  check  whether  it  contains  a  unique  subsegment  e' 
connecting  a  point  p^  on  8^  to  a  point  p2  on  82,  such  that  the  interior  of 
e'  is  fully  contained  in  R.  K  not  all  edges  in  S  have  such  vmique 
subsegments,  R  is  screened  out.  Otherwise,  sort  S  by  the  order  of  the 
endpoints  of  the  subsegments  e'  either  along  81  or  along  82,  to  obtain  a 
sequence  ^  of  edges. 

(IV)Output  the  collection  C{A,B)  of  all  sequences  ^  calculated  by  the 
preceding  step. 

It  is  easy  to  check  that  C(A,B)  contains  at  most  n  sequences,  because  at 
most  n  regions  R  can  contain  A  or  5  on  their  boundaries  (consider  the  fan  of 
the  n  curves  ^(a:)  emanating,  say,  from  A;  each  region  R  containing  A  on  its 
boundary  must  contain  one  of  the  sectors  between  two  curves  y(X),  yiX), 


adjacent  to  each  other  in  their  circular  order  around  A). 

Also,  the  preceding  considerations  clearly  imply  that  ^iA,B)  belongs  to 
C{A,B)  (this  is  easily  seen  to  be  true  even  in  case  5  is  a  ridge  point  relative  to 
A,  if  we  take  C(A,5)  to  be  any  of  the  edge  sequences  crossed  by  shortest 
paths  from  A  to  B). 

Our  next  goal  is  to  show  that  the  output  of  the  procedure 
EDGE_SEQUENCES  does  not  change  as  A  and  5  vary  slightly,  unless  one  of 
the  following  critical  conditions  occur: 

Definition:  A  critical  point  (A,B)  6  Z  is  defined  as  a  point  in  which  one  of 
the  following  critical  situations  occurs: 

(a)  Either  A  or  5  is  a  vertex  of  A"  or  is  a  ridge  point  relative  to  some  vertex 
ofK. 

(b)  Either  A  lies  on  'n(X,B)  or  B  lies  on  ir(X,A)  for  some  vertex  X  of  K. 

Proposition  2.4:  Suppose  (A, 5)  6  Z  is  a  non-critical  point.  Then  the 
procedure  EDGE.SEQUENCES  will  yield  the  same  output 
C{A'  ,B')  =  C(A^)  when  applK-4  to  any  pair  (A' ,5')  in  a  sufficiently  small 
neighborhood  of  (A,B). 

Proof:  First  note  that  the  conditions  of  criticality  define  a  closed  subset  of  Z, 
so  that  the  set  of  non-critical  points  is  open.  Consider  the  collection  of  curves 
y(X)  =  ^x.flW,  for  X  a  vertex  of  K.  We  first  claim  that  as  (A'  ,B')  varies  in 
a  sufficiently  small  neighborhood  of  (A,B)  which  consists  only  of  non-critical 
points,  the  curves  'Ya'.s'W  remain  simple  arcs  which  vary  continuously  with 
{A'  yB'),  their  points  of  intersection  vary  continuously,  and  the  combinatorial 
pattern  of  their  intersections  (that  is,  the  curves  intersecting  each  ^(X)  and 
the  order  of  their  intersections  along  y{X))  remain  constant. 

To  see  this,  note  first  that  as  long  as  neither  A'  nor  B'  becomes  a  ridge 
point  relative  to  some  vertex  of  K,  the  paths  tt(A',Z),  tt{B' ^)  vary 
continuously  with  A' ,  B' ,  for  each  vertex  X  of  K.  Furthermore,  their 
concatenation  is  a  simple  arc  unless  either  A'  lies  on  'n{B'  ^)  or  B'  lies  on 
'ir(A'^),  a  case  which  does  not  occur  as  long  as  (A'  JB')  remains  noncritical. 
U  two  arcs  yiX)  and  ^(7)  intersect  one  another,  say  'n{A' yX)  intersects 
iT(S',y),  then  since  (A',B')  is  assumed  to  be  noncritical  it  follows  that  these 
paths  must  intersect  one  another  transversally  so  that  their  intersection  point 
must  also  vary  continuously  with  A' ,  B' .  Finally,  note  that,  in  general,  if  a 
collection  of  simple  continuous  arcs  vary  continuously  in  such  a  way  that 
these  arcs  always  intersect  one  another  only  transversally,  then  the 
combinatorial  pattern  of  their  intersections  can  change  only  when  three  of 
these  arcs  intersect  at  a  common  point.  This  however  is  impossible  in  our 
case.  Indeed,  suppose  that  7a',5'(^),  7A',fl'(^)»  1a',b'{^)  intersect  at  a 
common  point,  for  some  three  distinct  vertices  X,  Y,  Z  of  K.  Then  two  of 
these  curves  must  intersect  at  their  portions  between  their  corresponding 
vertices  and  the  same  point  A'  or  B' ,  e.g.  'n{A' ,X)  intersects  7T(A',y).  But 
this  is  impossible  by  Lemma  2.1. 


-7- 

This  implies  that  the  combinatorial  structure  of  the  comiected  regions 
into  which  S  is  partitioned  by  the  arcs  7a',b'W  remains  constant  as  (A'  J3') 
varies  within  a  sufficiently  small  neighborhood  of  {A,B),  and  that  each  of 
these  regions  vary  continuously  with  (A'  JS'). 

Thus  the  first  step  of  the  procedure  EDGE_SEQUENCES  yields  the 
same  output  for  all  {A'  ,B')  sufficiently  near  (A^).  The  second  step  of  that 
procedure  calculates  the  sequence  of  edges  of  K  intersecting  each  half- 
boundary  5  of  regions  R  containing  both  A  and  B.  The  output  of  this  step 
can  in  fact  change  at  points  A' ,B'  for  which  two  curves  7A',fl'(-X'),  7a',s'(^) 
intersect  at  a  point  lying  on  an  edge  of.  K,  so  that  this  intersection  point  can 
newly  appear  in,  or  disappear  from  such  a  half -boundary  5 .  Nevertheless  we 
claim  that  these  potential  changes  cannot  affect  the  final  output  of  our 
procedure.  In  fact  we  need  to  show  that  for  each  region  R  =  Ra^b  in  the 
above  partitioning  of  S,  whose  boundary  contains  both  A  and  B,  so  that  these 
points  partition  it  into  two  half  boundaries  S^  =  8i(A,5)  and  82  =  82(A3), 
the  following  properties  hold: 

(a)  An  edge  e  intersects  both  8i(A,B)  and  82(A,5)  if  and  only  if  it  intersects 
both  8i(A'  ^')  and  82(A'  ,B')  for  all  points  (A'  ,B')  in  a  sufficiently  small 
neighborhood  of  (A,B). 

(b)  For  edges  e  as  in  (a),  there  exists  a  unique  subsegment  of  e  whose 
endpoints  pi  =  pi(AyB),  p2  =  p2iA,B)  lie  respectively  on  8i(A,j5), 
82(A^)  and  its  interior  is  contained  in  R,  if  and  only  if  this  property 
holds  for  all  (A' ^')  sufficiently  near  (A,5). 

(c)  For  edges  e  as  in  (a)  and  (b),  each  of  the  endpoints  pi{A'yB')  and 
PiiA'  ,B')  varies  continuously  with  (A'  ,B')  near  {A,B). 

Indeed,  suppose  for  the  moment  that  (a),  (b),  (c)  have  been  established. 
Then,  for  each  such  region  R,  asA',B'  vary  slightly  near  A,  5,  it  is  clear  that 
the  third  step  of  procedure  EDGE_SEQUENCES  will  produce  the  same 
sequence  of  edges  intersecting  both  8^  and  82,  that  if  this  step  does  not  screen 
out  R  at  (AyB)  it  will  not  do  so  at  nearby  points  (A'  ,B')  and  vice  versa,  and 
finally  that  the  points  pi(A'  JB')  and  piiA'  ,B')  calculated  by  this  step  will 
vary  continuously  along  hi(A' JB')  and  82(A',5')  respectively,  so  that  their 
order  along  these  half-boundaries  will  not  change.  Hence  the  collection  of 
edge  sequences  output  by  the  procedure  will  remain  constant  as  (A', 5') 
varies  sufficiently  near  (A ,5). 

To  prove  (a),  let  e  be  an  edge  of  K  which  intersects  both  portions  8^  and 
82  of  the  boundary  of  R,  and  let  e'  be  a  subsegment  of  e  whose  endpoints 
Pi(A,B),  p2{AJS)  are  intersections  of  e  with  81,  82  respectively,  and  whose 
interior  is  fully  contained  in  /?.  If,  say,  at  p-^{A,B)  the  edge  e  meets  just  one 
of  the  curves  7(X),  then  it  is  clear  that  as  A' ,  B'  vary  slightly  near  A,  5,  the 
point  pi(A',5')  will  continue  to  lie  on  ^{X)  and  will  vary  continuously 
(because  the  intersection  between  ^(X)  and  e  is  transversal).  Even  if  pi(A,5) 
lies  at  an  intersection  point  of  two  curves  7 (AT),  yiY),  it  is  still  easy  to  see, 
using  the  fact  that,  when  properly  unfolded,  y{X),  ^(K)  (and  e)  are  all  locally 
straight  segments  crossing  one  another  transversally  near  that  intersection 
point,    that    the   point  pi(A',5')    will    continue   to    exist,    and   will   vary 
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continuously  in  this  case  too  (possibly  moving  from  y{X)  into  7(7)  or  vice 
versa  as  A',  B'  vary  near  A,  B).  This  implies  in  particular  that  e  will 
continue  to  intersect  both  81  and  82,  for  points  A'  ,B'  near  A,B,  and  this 
clearly  implies  (a). 

To  show  (b),  note  that  the  condition  stated  there  can  change  only  when 
one  of  the  half-boundaries  81,  82  touches  e  without  crossing  it.  This  however 
is  impossible  because  the  intersection  of  any  -/(X)  with  e  is  always 
transversal,  unless  it  occurs  at  an  endpoint  of  e,  in  which  i:ase  it  cannot  cause 
changes  in  the  condition  stated  in  (b). 

Finally,  property  (c)  follows  directly  from  the  proof  of  (a)  given  above. 
The  proof  of  the  Proposition  is  thus  complete.  □ 

Proposition  2.4  implies  that  if  we  partition  the  subset  of  all  non-critical 
points  in  Z  into  connected  components,  then  within  each  such  component  the 
set  C{AyB)  of  the  0(n)  possible  values  that  i{A,B)  can  assume,  as  computed 
by  procedure  EDGE_SEQUENCES,  remains  constant.  To  estimate  the 
number  of  connected  non-critical  regions  in  Z,  we  show  next  that  the  critical 
points  in  Z  lie  on  the  union  of  O(n^)  critical  curves,  each  of  which  is  either  a 
straight  segment  or  part  of  an  hyperbola.  Hence  the  total  number  of 
intersections  between  these  critical  curves,  and  thus  the  maximal  number  of 
non-critical  regions  into  which  they  decompose  Z,  is  at  most  Oin'^),  so  that 
the  maximal  number  of  shortest-path  edge  sequences  induced  by  points 
(A  ,5)  e  ^1X^2  is  at  most  0(n^). 

To  see  this,  we  examine  each  condition  of  criticality  separately. 

Critical  points  of  type  (a).  A  critical  point  (A^)  of  type  (a)  is  one  for 
which  either  A  or  fi  is  an  endpoint  of  one  of  the  segments  ei,  ^2*  or  is  a  ridge 
point  with  respect  to  a  vertex  of  K.  As  mentioned  above,  each  vertex  of  K 
has  at  most  0{n)  ridge  points  on  ei  and  on  62-  Hence  critical  points  of  type 
(a)  lie  on  O(n^)  straight  segments  of  the  form  A  =  const  or  5  =  const  in  the 
parametric  space  Z. 

Critical  points  of  type  (b).  Fix  some  vertex  X  of  K,  and  partition  the  edge 
ei  into  0{n)  disjoint  subintervals  by  the  ridge  points  relative  to  X  lying  on  it. 
For  each  such  interval  /  the  sequences  ^(X,A),  for  A  6  /,  are  all  equal. 
Suppose  ^2  belongs  to  this  common  sequence,  and  unfold  the  fixed  sequence 
of  faces  of  K  crossed  by  all  paths  'Tr(X,A),  A  6  /,  such  that  all  these  faces  lie 
in  a  common  plane  and  such  that  after  this  unfolding  the  paths  7t(X,A) 
become  all  straight  segments.  Then  the  locus  of  pairs  (A ,5)  ^  Z  for  which 
A  €  /  and  B  lies  on  'tt(X,A)  is  equivalent  to  the  locus  of  all  pairs  (A* ,5*)  of 
points  lying  respectively  in  the  two  coplanar  unfolded  images  of  e^,  e2,  such 
that  the  line  connecting  A*  to  5*  passes  through  a  fixed  point  (the  unfolded 
image  of  X) .  An  easy  exercise  in  analytical  geometry  shows  that  this  locus  is 
a  hyperbola,  which  degenerates  into  a  straight  line  if  the  unfolded  images  of 
ei  and  ^2  are  parallel. 

Repeating  this  argument  for  each  vertex  X  of  K,  and  also  reversing  the 
roles  of  ei  and  ei,  we  conclude  that  type  (b)  critical  points  lie  on  O(n^) 
hyperbolae  and  straight  segments  in  Z. 
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This  analysis  establishes  the  claim  made  above. 

Finally,  note  that  to  estimate  the  total  number  of  shortest-path  edge 
sequences  along  S  it  suffices  (by  trimming  an  initial  and  a  final  portions  of 
shortest  paths  if  necessary)  to  consider  only  the  case  in  which  the  source  and 
target  points  of  the  shortest  paths  lie  on  edges  of  K.  Thus  by  repeating  the 
above  analysis  for  all  0(n^)  pairs  of  edges  of  K,  we  obtain  the  following 
main  theorem  of  this  section: 

Theorem  2.5:  The  total  number  of  shortest-path  edge  sequences  along  the 
surface  of  a  3-D  convex  polyhedron  with  n  vertices  is  at  most  0{nP). 

Remark:  We  do  not  know  whether  this  bound  is  tight.  In  fact,  the  best  lower 
bound  we  know  of  for  the  number  of  shortest-path  edge  sequences  along  S  is 
only  Cl{n}),  which  can  be  realized  even  for  a  fixed  source  point.  Our  intuition 
leads  us  to  conjecture  that  our  upper  bound  can  be  substantially  improved, 
and  an  obvious  open  subproblem  is  to  determine  whether  the  correct  bound  is 
indeed  only  O(n^),  or  else  to  find  an  example  with  more  than  quadratically 
many  shortest-path  edge  sequences. 

The  above  analysis  also  makes  it  easy  to  calculate  all  shortest-path  edge 
sequences  along  S,  in  the  following  crude  and  straightforward  manner: 

For  each  pair  e^,  ^2  of  edges  of  K  do  the  following: 

(a)  Draw  all  critical  lines  and  hyperbolae  in  Z  =  ^iXej  to  obtain  a 
decomposition  of  Z  into  connected  noncritical  regions. 

(b)  For  each  such  region  R,  choose  an  arbitrary  pair  {AJB)  €  R,  and 
apply  the  procedure  EDGE.SEQUENCES  to  (A  ,5),  obtaining  a 
collection  C{R)  =  C{A,B)  of  0{n)  edge  sequences. 

Output  the  union  of  all  collections  obtained  in  step  (b) . 

The  nmning  time  of  this  procedure  is  clearly  polynomial  in  /i.  In  fact,  it 
can  be  shown  that  step  (a)  can  be  accomplished  in  time  C>(n '^log  n)  per  each 
pair  ei,e2  of  edges,  and  step  (b)  can  be  implemented  in  time  C>(n^log  n)  per 
each  noncritical  region.  Hence  the  above  procedure  can  be  carried  out  in 
overall  C>(n ^log  n)  time. 

3.   Shortest  paths  amidst  a  fixed  number  of  convex  polyhedral  obstacles 

As  outlined  in  the  introduction.  Theorem  2.5  can  be  applied  to  obtain 
algorithms  for  the  calculation  of  shortest  paths  amidst  a  fixed  number  k  of 
convex  polyhedral  obstacles  having  altogether  n  vertices,  which  run  in  time 
polynomial  in  n  (albeit  exponential  in  k). 

Let  Ai,  .  .  .  ,Ajt  be  k  disjoint  convex  polyhedra  in  3-space,  and  let  X  and 
Y  be  two  points  in  free  space  (i.e.  in  the  complement  of  the  union  of  these 
polyhedra).  The  problem  is  to  calculate  the  shortest  path  tt  from  X  io  Y 
which  does  not  penetrate  into  the  interiors  of  A i,  ...  .A^.  As  shown  in  [SS], 
TT  is  a  polygonal  path  whose  comers  lie  on  edges  of  the  polyhedra  A/ ,  such 
that  with  each  edge  it  crosses  it  subtends  equal  incoming  and  outgoing  angles. 
Furthermore,  the  fact  that  each  A,  is  convex  implies  that  if  W,  Z  are  two 
points  on  it  which  lie  on  the  surface  of  the  same  polyhedron  A/  then  the 
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entire  subpath  of  it  between  W  and  Z  must  lie  along  that  surface.  Also,  each 
connected  portion  of  it  in  which  it  does  not  touch  any  polyhedron  must  be  a 
straight  segment  connecting  a  "take-off"  point  on  the  surface  of  one 
polyhedron  to  a  "landing"  point  on  the  surface  of  another  (with  X  as  an  initial 
take-off  point  and  y  as  a  final  landing  point).  It  is  also  easy  to  show  that  each 
of  those  take-off  and  landing  points  lie  on  a  polyhedron  edge  (unless  it 
coincides  with  one  of  the  points  X,  Y).  See  Fig.  5  for  an  illustration  of  this 
structure. 

Summing  up  these  observations,  there  exist  a  subsequence 
Ai^yA,^,  .  .  .  ,At  of  the  given  polyhedra,  with  each  polyhedron  A,  appearing  at 
most  once  in  this  sequence,  and  a  corresponding  sequence  of  take-off  and 
landing  points  li,ti,l2^h>  •  •  •  .'p»^p  with  /,,/^  lying  on  the  surface  of  A,^  for 
q  =  1,  .  .  .  ,p,  such  that 

IT  =  X/i  II  TTi  II   ri/2  II   TT2  II   ^I's  H     '   '   '     W   tp -lip    W   TTp    W   tpY 

where  each  tt^  is  a  shortest  path  along  the  surface  of  A,^  between  /,  and  t^ , 
and  where  the  interior  of  each  of  the  straight  segments  f,/,  +  i  does  not 
intersect  any  of  the  polyhedra. 

Hence,  to  find  the  desired  path  ir,  one  needs  to  solve  two  subproblems: 

I.      Determine  the  sequence  At^^i^,  .  .  .  ,A,^  of  polyhedra  through  which  ir 

passes. 

n.    Given  that  sequence,  calculate  the  take-off  and  landing  points  of  tt  on 
those  polyhedra. 

As  for  the  first  subproblem,  we  do  not  propose  any  efficient  procedure, 
and  make  do  with  enimieration  of  all  0(k^)  possible  subsequences  of 
polyhedra,  choosing  the  one  yielding  the  shortest  of  all  the  paths  calculated 
for  each  of  these  sequences.  (We  note  that  this  subproblem  is  at  least  as 
difficult  as  the  general  problem  of  calculating  shortest  paths  amidst  arbitrary 
polyhedral  obstacles,  for  which  only  exponential-time  algorithms  are  known 
as  yet  [SS],  [RS].) 

Suppose  then  that  we  are  given  a  specific  sequence  of  polyhedra,  which, 
without  loss  of  generality,  we  will  take  to  be  the  sequence  Ai^2.  •  •  •  --^t 
itself,  and  consider  the  second  subproblem  of  determining  the  take-off  and 
landing  points  on  the  shortest  path  from  X  to  y  which  is  constrained  to  pass 
through  the  surfaces  of  A^,  .  .  .  .A^  in  that  order.  Write  tt  as  above  as 
X/i  II  TTi  II  t]^l2  II  172  II  •  •  •  W  tk-ilk  W  TTi^  W  tj^Y,  whcrc  each  tt,  is  a  shortest 
path  along  the  surface  of  A,  from  /,  to  r,.  Note  that  by  Theorem  2,5,  for  each 
i  =  1,  .  .  .  ,k,  the  total  number  of  shortest-path  edge  sequences  through 
which  IT/  can  pass  is  at  most  0(n^),  where  n  is  the  maximal  number  of 
vertices  on  each  of  the  polyhedra  A,.  Hence  altogether  there  are  at  most 
0(n^*)  possible  sequences  of  edges  through  which  -it  can  pass. 

Let  4  be  one  of  these  candidate  sequences  of  edges.  We  first  unfold,  for 
each  I  =  1,  .  .  .  ,fc,  all  faces  of  A/  bounded  by  the  edges  in  ^  which  lie  on  A/, 
so  as  to  make  all  these  faces  co-planar,  and  to  make  tt/  a  straight  segment. 
These  unfoldings  then  reduce  the  specific  subproblem  at  hand  to  that  of 
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calculating  the  shortest  path  between  two  given  points  which  is  constrained  to 
pass  through  2k  lines  (namely  the  lines  containing  the  edges  in  ^  on  which  the 
take-off  and  landing  points  lie).  This  problem  however  is  solvable  in  time 
singly-exponential  in  k  [RS],  where  this  solution  can  be  obtained  using 
precise  symbolic  calculations,  allowing  one  e.g.  to  compare  in  precise  terms 
the  length  of  tt  with  the  length  of  any  other  shortest  path  constrained  to  pass 
through  another  sequence  of  polyhedra.  (After  calculating  each  such  tt,  we 
also  have  to  check  that  the  take-off  and  landing  points  do  indeed  lie  within 
their  corresponding  edges,  and  not  just  on  the  lines  containing  those  edges.  If 
this  does  not  hold,  tt  is  discarded  as  a  candidate  shortest  path.) 

Finally,  note  that  for  each  candidate  path  tt  we  also  need  to  check  that  it 
does  not  intersect  any  polyhedron  other  than  those  through  which  it  is 
assumed  to  pass  a  priori,  and  consider  only  paths  tt  which  meet  this  criterion. 
This  additional  test,  however,  is  easily  accomplished  in  polynomial  time,  e.g. 
by  testing  each  segment  of  tt  for  intersection  with  each  of  the  remaining 
polyhedra  faces. 

We  thus  obtain  the  following  main  theorem  of  this  section: 

Theorem  3.1:  The  shortest  path  between  two  points  X,  Y  which  does  not 
penetrate  into  the  interiors  of  k  disjoint  convex  polyhedra  having  altogether  n 
vertices,  can  be  calculated  in  time  n'^^*). 

Remark:  This  result  is  related  to  the  solution  of  the  2-D  "Asteroid  avoidance 
problem"  studied  in  [RSh].  Theorem  3.1  is  however  stronger  than  that  result, 
because  it  makes  use  of  the  rather  heavy  Theorem  2.5.  In  contrast,  the 
analysis  in  [RSh]  deals  with  only  trivial  sets  of  edge  sequences,  whose 
calculation  is  also  straightforward. 
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